A Passive Phase Noise Cancellation Element 



O 

(N 



m 



Eyal Kenig,!' * M. C. Cross/ Ron Lifshitz,^ R. B. Karabalin/ L. G. Villanueva/ M. H. Matheny/ and M. L. Roukes^ 

^ Kavli Nanoscience Institute and Condensed Matter Physics, 

California Institute of Technology, MC 149-33, Pasadena, California 91125, USA 

^Raymond and Beverly Sackler School of Physics and Astronomy, Tel Aviv University, 69978 Tel Aviv, Israel 

(Dated: March 14, 2012) 

We introduce a new method for reducing phase noise in oscillators, thereby improving their 
frequency precision. The noise reduction device consists of a pair of coupled nonlinear resonating 
elements that are driven parametrically by the output of a conventional oscillator at a frequency 
close to the sum of the linear mode frequencies. Above the threshold for parametric response, the 
coupled resonators exhibit self-oscillation at an inherent frequency. We find operating points of the 
device for which this periodic signal is immune to frequency noise in the driving oscillator, providing 
a way to clean its phase noise. We present results for the effect of thermal noise to advance a broader 
understanding of the overall noise sensitivity and the fundamental operating limits. 

PACS numbers: 05.45.-a, 84.30.-r, 85.85.+J, 62.25.-g. 
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The emergence of self-oscillation has a major scien- 
tific significance as a widespread phenomenon in physics, 
chemistry, and biology [1] . Oscillators are also extremely 
useful, frequently appearing as crucial elements in the 
electrical devices that surround us in our highly techno- 
logical environment. Essentially, oscillators are devices 
generating a periodic signal at an inherent frequency, 
whose primary function is therefore to provide a time or 
a frequency reference. An ideal self-sustained oscillator 
is mathematically described as a limit-cycle in the phase 
space of dynamical variables, or equivalently as a periodic 
solution of a set of autonomous differential equations, in- 
dependent of an external time reference. The ideal os- 
cillator can thus be described in terms of a steadily in- 
creasing phase variable corresponding to the phase space 
point advancing around the limit-cycle, with a 27r phase 
change corresponding to a period of the motion. This 
phase is highly sensitive to additional stochastic terms, 
or noise, in the equations of motion, as the appearance of 
periodicity without an external time reference implies the 
freedom to drift along the phase direction. The stochas- 
tic phase dynamics lead to a broadening of the peaks in 
the power spectrum of the oscillator output, which arc 
perfectly discrete in the ideal case [2, 3], and to a degra- 
dation of the performance as a time or frequency refer- 
ence [4]. Thus, an essential task in the design of a good 
oscillator is to reduce the effects of the noise, present in 
the system, on the oscillator phase. 

In this letter we propose a general scheme to reduce, or 
even eliminate, the noise in the output of an oscillator by 
passing the signal through a second passive noise cancel- 
lation device, rather than manipulating the oscillator it- 
self. It is therefore broadly applicable to enhance the per- 
formance of existing oscillator designs. It is also valuable 
from a basic physics perspective, eliminating the need to 
analyze the active resonator-amplifier feedback system: 
instead, the noise performance is mapped to a passive el- 
ement, whose fundamental stochastic properties are more 



amenable to the powerful tools of statistical physics [5] . 
The context of our work is to use nanoscale or microscale 
resonators to build high precision oscillators as illustrated 
in Fig. 1, although our scheme applies more generally. 
The nano- or micro-electromechanical systems (NEMS or 
MEMS) implementation shown in Fig. 1 consists of a pair 
of coupled resonating elements that are parametrically- 
driven with a noisy frequency near the sum of their lin- 
ear mode frequencies. As explained below, it produces 
a signal with reduced noise at a frequency near the dif- 
ference of these linear frequencies. The generation of a 
two-frequency signal by parametric excitation at the fre- 
quency sum is called non-degenerate parametric excita- 
tion, a phenomenon that has been studied in mechanical 
systems [6] , but is more common in the context of optical 
parametric oscillators [7], where our scheme also applies. 

The ability of an oscillator based on a nonlinear res- 
onator to suppress phase noise was demonstrated by 
Yurke et al. [8] , who studied a nonlinear beam resonator 
driven into self-oscillation through an active feedback 
loop, composed of an amplifier driven to saturation and 
a phase shifter. They showed that if the resonator is 
operated at the critical Duffing point, by applying the 
minimal drive required for making the amplitude vs. fre- 
quency curve multivalued, and by choosing the phase 
shift in the feedback loop between the signal and the 
drive to be 27r/3, the noise in the phase of the fedback 
signal has no effect on the oscillator phase. This was 
understood in terms of the insensitivity of the resonator 
frequency to the drive phase at this specific operating 
point. Here we generalize this understanding of noise re- 
duction in the feedback oscillator, and show that in the 
present case the output phase noise due to phase noise in 
the drive may be eliminated by operating at any points 
for which the output frequency itself is insensitive to the 
drive frequency. 

We initiate our analysis by modeling a pair of coupled 
resonators, like the ones shown in Fig. 1, by two dinien- 
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FIG. 1. (Color online) An illustration of the phase noise can- 
cellation scheme. An oscillator produces a signal with a noisy 
frequency around oii -|-a;2. This signal parametrically drives a 
pair of coupled NEMS or MEMS beams with a relative phase 
of 180° . The output signal at a frequency around aj2 — ^i is 
given by squaring and filtering. 
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where wj = 1 and a;2 = 1 + 2A arc the linear frequen- 
cies of the in-phase and out-of-phase modes, respectively, 
eflp ~ LOp — uji — 102 is the small difference between the 
drive frequency and the sum of the linear mode frequen- 
cies, and g ~ {hi — /i2)/2. Similar slow equations for 
pairs of coupled resonators were recently introduced to 
study chaotic dynamics [11] and to analyze the so-called 
Bifurcation- Topology Amplifier [12]. 



In magnitude-phase coordinates. An 
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Eqs. (2) can be transformed into four dynamical equa- 
tions for the variables ai, 02, $ = 01+02, and ^I* = 



sionless equations of motion 



in +Xn + xl + tin + r\x\xn + thnXn COs{u;pt) 



-1- A(a;„ - Xk) = 0, 



(1) 



with (ri, fc) ~ (1,2) and (2,1). For details on the use of 
such equations for modeling NEMS and MEMS devices 
see Lifshitz and Cross [9, 10]. We only wish to highlight 
the following points. The resonators are taken to be iden- 
tical, having a resonant frequency wq, which has been 
scaled out, although considering different resonant fre- 
quencies would not qualitatively change our results. The 
nonlinearity of the resonators originates from both the 
elastic restoring force and damping mechanism. In micro- 
and nano-scale resonators the linear damping is typically 
weak, e ^ 1, where 1/e is the quality factor of the res- 
onator, and correspondingly a small drive amplitude is 
sufficient to excite them. Accordingly, each resonator is 
parametrically excited with a drive amplitude e/i„. In 
this regime, we can focus on the slow-time modulation 
of the basic oscillatory motion of the in-phase and out- 
of-phase modes of the coupled resonators, described by 
a pair of complex equations for the corresponding mode 
amplitudes An{T), where T = ef is a slow time scale. 
In the regime of strong coupling, where the linear mode 
splitting is much larger than the resonator bandwidth. 
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The right hand sides of these equations depend only on 
the three variables oi , 02 and $. The absence of a depen- 
dence on \1/ reflects the fact that Eqs. (2) are unchanged 
by the transformation (01, 02) — >■ (0i + /?, 02 — /?)■ This 
property of the non-degenerate response to parametric 
excitation is well known [13], and does not occur in the 
degenerate case, where wi = W2, in which both phases 
are fixed, corresponding to oscillations that are locked 
to the phase of the drive. The fixed point solutions 
(oi.o, 02,0, ^0) of the three dynamical equations (3a)-(3c) 
correspond to periodic orbits of Eqs. (2) with frequency 
rio/2, where 0.q = /4(ai_o, 02,0, *&o)- The amplitudes and 
frequency of these solutions are shown in Fig. 2 as a func- 
tion of the drive frequency for g = 10, 77 = 1, and a cou- 
pling A = 7, corresponding to a linear mode splitting of 
\/l + 2A — 1 ~ 2.87 in units of the resonance frequency 
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FIG. 2. The periodic solutions of Eqs. (2). (a) The squared 
mode amplitudes (the in-phase mode has the larger ampli- 
tude); (b) Twice the frequency of the periodic solutions. Solid 
and dashed lines are stable and unstable solutions, respec- 
tively. The parameters are A = 7, g = 10, and r; = 1. 
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We now demonstrate how the dynamics can be utilized 
to eliminate the phase noise of the driving oscillator. We 
model this noise as white noise of intensity Fq in the 
frequency of the drive, replacing ftp with flp + 5p(T), 
and assuming that {Ep{T)Ep{T')) = Fn^S(T-T'). This 
noise causes the phase of the driving source $p to diffuse 
as {[<i>p{T + r) - ^p{T)]^) = -FoJt|. Additional noise in 
the amplitude of the drive signal can be suppressed by 
using a limitcr or a highly saturated amplifier [8, 14]. The 
output signal, which is obtained by mixing and filtering 
the response to this drive, is proportional to cos[(w2 — 
uJi)t + ^{T)]. At the steady state of oscillation, the phase 
^ can be expressed as ^(T) = floT+ifj{T), where ip is the 
small stochastic perturbation induced by the noise Sp. 
The phase perturbation ijj can be calculated by solving 
the linearized version of the polar amplitude equations 
(3) spectrally, as was done for a single oscillator by Yurke 
et al. [8], or using the secular perturbation method of 
Demir et al. [3] . This gives diffusion of the phase of the 
output signal, with the variance of the phase difference 
growing linearly in time 



mT + T)-i;{T)f)^Dn,FnjT\, 



(4) 



where the diffusion constant Dq quantifies the sensitiv- 
ity of the phase of the output signal to the noise Ep in 
the drive frequency. This expression applies in the limit 
that the time r is much longer than the decay time onto 
the limit cycle for the noiseless oscillator. This approxi- 
mation corresponds to frequency offsets close to the os- 
cillator frequency. This is the relevant regime to explore 
since the width of the oscillator spectral peak is typically 
narrower than that of the driven resonator, whose width 
is determined by the quality factor, which in turn sets 
the decay time. The phase diffusion induced by white 
drive-frequency noise corresponds to a Lorenzian spec- 
tral peak of the output signal with a width of Dq Fq 
times the resonator bandwidth eujQ. Equation (4) is the 
expected result because after transients have decayed and 
the system has settled onto the limit cycle, the net effect 



FIG. 3. Phase diffusion of the output signal induced by fre- 
quency noise of unit strength for different values 77 of the 
nonlinear damping. For each curve, SI™'''' is the value of Qp 
for which the diffusion is zero. Note that Dq^ is also smaller 
than the phase diffusion of the driving oscillator along most 
of the curve, and that the solid curve is the squared derivative 
of the solid curve in Fig. 2(b). 



of the stochastic dynamics is a Brownian motion of the 
free phase. 

Without going into the details of the calculation, we 
note that the diffusion coefficient D^ for a white noise 
source N in Eqs. (3) may be calculated as the squared 
scalar product of two vectors Dat = (vj_ • vjv)^. The 
vector vj_ captures the phase sensitivity of the system, 
through the Jacobian matrix of Eqs. (3) describing the 
linearized flow in the vicinity the limit-cycle. Specifically, 
it is the eigenvector of the transpose of this Jacobian ma- 
trix that corresponds to the zero eigenvalue, vat is the 
noise vector, whose nth entry is the term multiplying the 
noise source in the function /„ (vw = (0,0, 1,0) for the 
current example since flp only appears in Eq. (3c) and 
it has a coefficient of magnitude one). This description 
of the diffusion coefficient is a simplification of the gen- 
eral approach described in Ref. [3] for systems in which 
motion along the limit-cycle is described by a steadily ad- 
vancing phase variable ^ which does not affect the time 
evolution of the dynamical system, and thus the limit- 
cycle is represented by a fixed point in all other variables. 
In terms of the current example, this reducibility of the 
dynamical system is expressed by the fact that ^ does 
not appear on the right hand sides of Eqs. (3). In the 
more general case the corresponding vectors vx(i) and 
V7v(i) are time dependant: they are periodic, having the 
period of the limit-cycle, and the diffusion coefficient is 
given by the time average of their squared scalar product. 

In systems such as the current one under consideration, 
calculating the zero mode of the transposed Jacobian ma- 



trix provides a way to calculate the phase diffusion that 
results from any white noise vector acting on a limit- 
cycle. However, if the noise originates from fluctuations 
in some parameter of the equations p so that the noise 
vector is VAT = Vp = (^, ^, •••, -^j, the long time 
phase diffusion is directly related to the dependence of 
the oscillation frequency on the parameter p through 
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The second equality follows from a perturbation analy- 
sis of the change in frequency 5^p due to a small per- 
turbation 5p in the parameter p. Eq. (5) shows that 
to reduce the frequency stability degradation due to pa- 
rameter noise we seek extremum points in the curves of 
the oscillation frequency vs. the noisy parameter. This 
possibility of complete noise elimination is due to the re- 
ducibility of the dynamical description, which makes the 
two vectors used for calculating the diffusion coefficient 
constant in time. In this case, it is possible to tune a 
single parameter and make these vectors orthogonal [15]. 

We now apply the ideas of the previous paragraphs to 
our noise reduction setup, for which the noisy parame- 
ter \s p — Vtp. Figure 2(b) shows an extremum we can 
exploit at the maximum of the frequency curve. By op- 
erating at this point, i.e. at the drive frequency fi™^'', 
the effect of drive frequency noise on the phase diffusion 
is eliminated, as shown in Fig. 3, and the corresponding 
frequency spectrum is narrowed (to a sharp peak in the 
absence of other noise terms). It is also noteworthy that 
for the parameters used to plot the diffusion coefficient 
Dq_^ in Fig. 3, the frequency stability is also improved 
along most of the curve, since Dq < 1 . and thus the out- 
put frequency noise is less than input frequency noise. An 
examination of the different curves in the figure indicates 
that larger nonlinear damping makes fl^^^ more acces- 
sible experimentally, since this increases the frequency 
separation between the maximum {dil.o/dnp = 0) and 
the saddle-node (dfJo/f^^p = oo). It is possible to con- 
trol the nonlinear parameters of a resonator (both Duff- 
ing and damping terms) by adding a parametric feedback 
loop, as we have recently shown elsewhere [16]. 

The ability of the device to clean phase noise is lim- 
ited by the thermal noise floor. To discover this fun- 
damental limit we add the complex thermal noise terms 
S„ = Sj^n +i'E.in to each of Eqs. (2), with the individual 
noise components white, uncorrelated, and of the same 
intensity Fth. Our analysis shows that the phase diffu- 
sion resulting from thermal noise is then given by Eq. (4), 
but with the coefficients DQ^Fa^ replaced by Dth-Fth, 
where Dth is plotted in Fig. 4. In terms of the actual 
physical parameters, the lower phase diffusion limit is 
(kBT/Ec)Dth (in units of the bandwidth ewo) with fceT 
the thermal energy, Ec = ern^ujQ/a the potential energy 
of the resonator element when driven to the Duffing crit- 




FIG. 4. Thermal noise limit of the phase diffusion of the 
output signal, for fth = 1 and the same parameters as in 
Fig. 2. The asterisk indicates the value at fi™^''. 



ical amplitude, a the Duffing nonlincarity parameter (in 
units of force per volume) , and m the resonator mass (see 
Ref. [10]). Since thermal noise does not originate from 
some fluctuating parameter, it cannot be eliminated by 
finding extremum points in the dependence of the out- 
put frequency. Nevertheless, it is possible to lower the 
thermal noise limit by reducing Z?th. The idea is based 
on the fact that the direct effect of thermal noise on the 
phase variable is reduced at large amplitudes, so in the 
large amplitude limit, the only cause of phase diffusion is 
the conversion of noise from the amplitude to the phase 
(AM-PM conversion). For a single oscillator, this conver- 
sion is completely eliminated by operating at a point for 
which the resonator frequency is insensitive to the ampli- 
tude [17], and this approach can be generalized to achieve 
amplitude-phase detachment in our system as well. 

In conclusion, we have described a passive device that 
eliminates phase noise in oscillators. The device is made 
of two coupled resonators, driven parametrically in the 
non-degenerate mode by the output of a conventional 
oscillator. We find a driving frequency for which the 
resulting limit-cycle oscillation frequency is insensitive to 
the drive frequency and show that operating at this point 
eliminates the phase noise in the driving oscillator. We 
have discussed the operational limitation due to thermal 
noise, and have suggested ways to improve this limit. 
Along with the interesting physics it portrays, this device 
offers a practical way to handle the extensively studied, 
cardinal problem of oscillator phase noise. 
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